Absorption and Emission in the non-Poisson case. 



O 

o 



Gerardo Aquino^, Luigi Palatella^, Paolo Grigolini^'^^'^ 

^Center for Nonlinear Science, University of North Texas, P.O. Box 311427, Denton, Texas 76203-1427 

^ Dipartimento di Fisica dell'Universita di Pisa and INFM, via Buonarroti 2, 56126 Pisa, Italy and 

^ Istituto dei Processi Chimico Fisici del CNR Area della Ricerca di Pisa, Via G. Moruzzi 1,56124 Pisa, Italy 

(Dated: February 2, 2008) 

This letter adresses the challenging problems posed to the Kubo- Anderson (KA) theory by the 
discovery of intermittent resonant fluorescence with a non-exponential distribution of waiting times. 
We show how to extend the KA theory from aged to aging systems, aging for a very extended time 
period or even forever, being a crucial consequence of non-Poisson statistics. 
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In the last few years, as a consequence of an increas- 
ingly faster technological advance, it has become clear 
that the conditions of ordinary statistical mechanics as- 
sumed by the lineshape theory of Kubo and Ander- 
son (KA) 1], are violated by some of the new mate- 
rials. For instance, the experimental research work of 
Neuhauser et al. |3] has established that the fluorescence 
emission of single nanocrystals exhibits interesting inter- 
mittent behavior, namely, a sequence of "light on" and 
"light off" states, departing from Poisson statistics. In 
fact, the waiting time distribution in both states is non- 
exponential, and it shows a universal power law behavior 
|3J. In this paper, for simplicity, we assign to both states 
the same waiting time distribution 
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with 1 < 1/ < cxD. The parameter T > is introduced for 
the purpose of making ^{t) finite at t = so as to ensure 
its normalization. We shall focus on the case when v <?,. 
In accordance with Brokmann et al. |j| , the experimental 
condition v < 2 implies that the observed waiting time 
distribution depends on the time at which observation 
begins. Let us assume that the probability of the first 
jump from the "on" ("off") to the "off" ("on") state is 
given by Eq. (Q, if the observation begins at i = 0. If the 
observation begins at a later time ta > 0, the distribution 
of the sojourn times, before the first jump, turns out to 
be different from Eq. (^: it is ta dependent and, for 
this reason, is denoted by f{t,ta). This is the property 
responsible for the breakdown of the ordinary KA theory: 
it is the aging effect on which we focus our attention in 
this letter. It is worth noticing that when v > 2, this 
aging effect is still present, in a less dramatic form, given 
the fact that a stationary condition exists, even if the 
regression to it takes a virtually infinite time if ;/ < 3 j^ . 
The authors of Ref. ^ showed how to derive the ab- 
sorption lineshape in the case ly > 2, when the stationary 
condition applies, and evaluated the form that the spec- 
trum would have, immediately after switching on the ra- 
diation field, when the non-stationary condition v < 2 
holds true. Here we illustrate a way to evaluate the time 



evolution of the absorption spectrum, so as to take into 
account the aging effects of Brokmann et at, with v < 2, 
as well as those of Ref. s^], with z^ > 2. We use the fol- 
lowing stochastic equation: 
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The quantity ^(t) is a complex number, corresponding 
to the operator |e)(ol of the more rigorous quantum me- 
chanical treatment |3, |e) and \g) being the excited and 
the groundstate, respectively, wq is the energy difference 
between the excited and the ground state, and ^(i) de- 
notes the energy fluctuations caused by the cooperative 
environment of this system. In the presence of the coher- 
ent excitation, Eq. |(2Jl becomes 



dt 



fi{t) = i{uJo + £,{t))fi{t) + kexp{iujt), 



(3) 



where lo denotes the radiation field frequency. It is 
convenient to adopt the rotating-wave approximation. 
Let us express Eq. by means of the transformation 
jl{t) = exp{iLut) ii{t) . After some algebra, we get a simple 
equation of motion for fl{t) . For simplicity we denote il{t) 
with the symbol fi{t) again, thereby making the resulting 
equation read: 



d , , 



z{5+at))Kt) 



(4) 



where 5 — luq ^ ^. The reader can easily establish the 
connection between this picture and the stochastic Bloch 
equation of Ref. 7] by setting fi — v + iu. Note that 
the three components of the Bloch vector in Ref. |3], 
(m, V, w), are related to the rotating-wave representation 
of the density matrix p, v and u being the imaginary 
and the real part of e~''^Vse, and w being defined by 
w = {pee — Pgg)l2-- Note that the equivalence with the 
picture of Ref. [7| is established by assuming the radia- 
tive lifetime of the excited state to be infinitely large and 
the Rabi frequency fi = fc vanishingly small. 

It is straightforward to integrate Eq. Q , thereby get- 
ting 



Pi{i)^k / dtV-f*'«(*")'**"e'*(*-*'), 



(5) 



with the dipole /i = 0, when the exciting radiation is 
turned on. Now we have to address the intriguing issue 
of averaging Eq. (|SJ) over a set of identical systeins,jii 
such a way as to take aging effects into account 
In fact the averaging process turns Eq. jSJ) into 
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with the subscript t' denoting that the system, brand new 
at i = 0, is t'-old when we evaluate the corresponding 
characteristic function, thereby implying that the dis- 
tribution of waiting times before the first jump, is not 
ipit), and f{t,t') has to be used instead. Our numeri- 
cal approach rests on the following prescriptions. Using, 
the distribution of Eq. |^ , we run N distinct sequences 
JTi}, with N ^ 1. For any of these N sequences the 
sojourn in one of the two states begins exactly a time 
t = and ends at time t = ti. For N/2 of these se- 
quences we use the "light on" as initial condition, and 
for N/2 the "light off' state. Let us consider the for- 
mer type of trajectories for illustration purpose. The 
"light on" state begins at t = and ends at i = ti, at 
which time the "light off' state begins, ending at time 
t = Ti + T2, and so on. With this numerical method 
we produce a set of fluctuations ^(i). Then, we create 
a set of diffusion trajectories J^, (,{t")dt" , and hence the 

set of exponentials e*-'*' '>'■* •'''* . Since all trajectories of 
this set are i'-old, the resulting numerical average is au- 
tomatically equivalent to evaluating {e^^t' ^(* ^'^* )t' with 
the waiting time distribution f{t,t'). We point out that 
the number of photons emitted at time t is determined 
by N{t) = (/i(i)/i*(t)}. It is straightforward to prove 
that the rate of photons emitted, namely, R{t) = dN/dt, 
obeys the relation 



R{t) ^ 2kRe{fi{t)). 



(7) 



Thus, we conclude that the real part of /x(i) can be used 
to denote either emission or absorption at time t. We 
note that the approximation ensuring the equivalence be- 
tween our picture and Ref. (7|, for large photon count, 
has also the effect of making the absorption identical to 
the emission spectrum. Fig. Q] illustrates the rate of ab- 
sorption changing upon change of time, when v = 2.5. 
We see that the spectrum changes from a Lorentzian 
shape centered at lo = loq to a bi-modal shape at later 
time, in accordance with earlier theoretical work |lflllll |. 
The authors of Ref. IQ] studied the aged condition 
f{t, oo) with 2 < v < 3, and found that the correspond- 
ing spectrum has two sharp peaks. Klafter and Zumofen 
proved that the same case with f{t,ta = 0) = ip{t) 
yields Levy diffusion, thereby implying an exponential 
characteristic function, and consequently a Lorentzian 
spectrum. It is remarkable that our theoretical perspec- 
tive establishes a connection between the prediction of 
Zumofen and Klafter, valid at short times, and that of 
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FIG. 1: The numerical (top) versus the theoretical (bottom) 
absorption spectrum with v = 2.5, T -W = 0.4. We assign to 
the fluctuation ^(f) the intensity W = 1. 



Ref. [l(|, valid at large times. Remarkably, we evalu- 
ate numerically also the spectrum time evolution in the 
case V < 2, where the process of aging keeps going on 
forever without ever reaching the stationary condition, 
which does not exist, in this case (see Fig. |5}. 

Let us show now how to reproduce these numerical re- 
sults with a proper theoretical treatment. Let us consider 
as an example a single trajectory starting a time t' with 
^ — -\-W and ending at time t with the same positive 
value after n switches of the variable ^ between its two 
values ±.W . In this case the integrand of (0, has the 
following form: 



„(a+(ti-t')+a-(*2-tl) + -|-a+(t-t„)) 



(8) 



where a± = i{5 ± W) and n can only be even. 

To determine the contribution to {ii{t)) of Eq. © 
stemming from this kind of trajectories, we have to av- 
erage the term ^ on the set of all possible sequences of 
the variable ^(i) running from t' to t with n switches, to 
sum over n and then to carry out the integration on t' . 
To this purpose, as earlier stated, besides the probabil- 
ity density ^^{t) of having an interval r between two 
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FIG. 2: The numerical (top) versus the theoretical (bottom) 
absorption spectrum with v = 1.5, T -W = 0.4. We assign to 
the fluctuation ^(i) the intensity W — 1. 



generic consecutive switching events, it is necessary to 
use f{T,t'), this being the conditional probability den- 
sity that, fixed a time t' , the first next switching event of 
the variable £,{t) occurs at time t' + t. The result of this 
procedure turns the contribution to @ into: 
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. pa+(t3-*2) . . . 
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where, in addition to the crucial probability density 
f{T,t'), we have used 



*(r) = 1 - / i){T)dT\ 
Jo 



(10) 
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FIG. 3: Comparison between numerical and theoretical ab- 
sorption spectra for diflerent values of the de-tuning param- 
eter 5 and of the power index v. Moving from the top to the 
bottom in the right hand portion of the figure, 5 = 1, v = 1.5, 
5 = 1, V = 2.5, 5 = 0.6, V = 2.5, 5 = 0.6, v = 1.5. The 
numerical curves are denoted by black circles, open circles, 
open triangles and black triangles, respectively, and in this 
scale are almost indistinguishable from the full lines denoting 
the corresponding theoretical results. 



which is the conventional probability that no switch oc- 
curs for a generic interval of time t, and 



F(T,i') = l- r f{r',t')d7 
Jq 



(11) 



which is the corresponding aging property, depending on 
f{T,t'), and indicating therefore the conditional proba- 
bility that, fixed t', no switch occurs between t = t' and 
t = t' + T. The overall factor of 1/2 of the contribution 
© is a consequence of the fact that at time t' the fluc- 
tuation ^(t), supposed to be positive, can get with the 
same probability the negative value. Let us address now 
the problem of finding an exact analytical expression for 
the crucial property /(t, <'). The exact expression for 

/(T,t'), is 



where 



/(r,0= / dr'G{t'-T')^{r + r'), (12) 
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It is straightforward to find the Laplace transform of 
G{t). This is given by 



CO ^ 
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(14) 



where V'(w) denotes the Laplace transform of ip{t). Thus, 
the Laplace transform of H12() with respect to t' reads: 



f{r,u') 



1 - V(w') 
thereby yielding: 







(15) 
(16) 



With this expression the prescription necessary to evalu- 
ate the contribution to (/i(i)) of Eq. ^ of the trajectories 
beginning in the "light on" state at t' and ending in the 
same state at t, is completed. 

As a last step, we calculate the Laplace transform of 
0, and of the equivalent expressions for all the other 
possible conditions of motion from t' to t, namely, with 
the noise ^(t) moving from +W to —W, from —W to +W 
and, finally, from —W to —W, as well. This procedure 
yields, as a final result, the Laplace transform of R{t), 
denoted by R{u), which is proved to have the following 
analytical expression: 

R{u) = kCMt) + ^^*m = ^{A+iu) + A_{u) + C{u)). 

(17) 
The explicit expression for Aj^{u), which is calculated 
taking into account the contribution to {^{t)) of those 
trajectories ending at time t with a positive value ^(t) — 
+W for the flucuating variable, is 

A^iu) = h F+[u), 



1 — ^p{u — a+)'ip{u — a_) 



(18) 



with: 



^>{u~a±)^C[^{t)e"^']^^-^^^^^-^, (19) 
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and 



F±{u) = C[ dt'F{t-t',t')e 



a±(t"t')l — 



l/u-f±{u) 

(21) 
A- (u) takes into account the contribution of all trajecto- 
ries ending at time t with the negative value ^(i) = —W, 
for the fluctuating variable and is derived from the ex- 
pression for A+{u) by replacing everywhere a± with aip 
(sending f±,F± ^ hJ^)- Thus, k{A+{u) + A.{u))/2 
represents the laplace transform of {fi{t)). Finally 
kC{u)/2 is the Laplace transform of (/i*(t)), and it is de- 
rived from the earlier expression for k{A^{u) + A^{u))/2 



by replacing everywhere a± with —a±. To establish a 
comparison with the result of the numerical experiment 
we have anti-Laplace transformed the analytical expres- 
sion of Eq. (|17l) using a Talbot algorithm implemented 
on Mathematica 5.0. The result of this procedure is il- 
lustrated by both Fig. Q]and Fig. |21 where the analytical 
predictions are compared to the corresponding numer- 
ical experiments. The overall qualitative agreement is 
remarkably good. Apparently, in Fig. ^the departure of 
the theoretical calculation from the fine structure of the 
numerical result is larger than in Fig. [3 Actually, the 
numerical algorithm produces fluctuations of the same 
intensity in both cases, although the more detailed scale 
of Fig. n makes them appear larger. To establish be- 
yond any doubt also the quantitative agreement, in Fig. 
|31 we compare the two predictions for specific values of 
the de-tuning parameter 6. This shows that also the 
quantitative agreement is very satisfactory. Furthermore, 
when a stationary condition exists, it is straightforward 
to prove that Eq. (|17|l recovers, as asymptotic-time limit, 
the same result as that of p . 

In summary, we could afford a complete picture of 
the aging process, going beyond either the short-time 
or asymptotic-time limitations of Ref. p. This letter 
shows also that these aging spectra are a consequence of 
the fact that the field acting at time t' makes the pho- 
ton emission depend on i'-old trajectories. This yields 
Eq. ©, with the subscript t' working as an age indica- 
tor. This formula intertwines the system to the radiation 
field, thereby violating the ordinary linear response pre- 
scription, which is recovered in fact in the Poisson case. 
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